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Abstract. The effect of the Zeeman field on the defect-induced Kondo effect in graphene is 
investigated. The effective model of the Kondo effect is derived, and is analyzed on the basis of 
the numerical renormalization group method. It is found that, under the Zeeman field, at the 
border of the spin polarized state and the Kondo- Yosida singlet state, there is an unusual fixed 
point where the entropy of the defect is Ub In 2 and the expectation value of Si is 1/8. 



1. Introduction 

The effects of defects on the properties of graphene are of scientific and practical importance pQ. 
For studying the effects of defects experimentally, ion irradiation has been used, which introduces 
point defects in graphene. Recently, in such an ion-irradiated graphene, the Kondo effect is 
observed [2]. Usually, the Kondo effect indicates the presence of localized moments in the 
sample. Therefore the above result suggests that a point defect in graphene causes a localized 
moment. The magnetic behavior of the defects in graphene is also observed in the recent 
magnetization measurement of ion-irradiated graphene [3] . 

As a possible origin of this Kondo effect, we proposed a localized moment on a sp 2 orbital 
around the defect [4j. There are three sp 2 orbitals around a point defect. With Jahn- Teller 
distortion, two of them form a covalent bond and remaining one behaves as a localized moment 
with spin-half. Furthermore, with this lattice distortion, there is a hybridization between the 
localized sp 2 orbital and conduction tt electrons, and the Kondo effect can occur. On this 
scenario, we introduced an effective model and analyzed it by the numerical renormalization 
group (NRG) method [6]. As a result, we showed that the experimentally observed gate 
voltage dependence of Kondo temperature can be understood by assuming that the hybridization 
between the sp 2 and tt orbitals is strong enough [JJ. 

In the above experiment, the negative magnetoresistance (MR), which is the evidence of the 
Kondo effect, was also observed. Although the characteristic magnetic field for this negative 
MR is found to be smaller by one order than the usual Kondo-effect case, its origin has not been 
understood. So far, the studies on the effect of magnetic field on the Kondo effect in graphene 
are limited, and there are few studies on the effect of Zeeman field [7]. 

Motivated by this negative MR, in this paper, we construct an effective model of graphene 
with a point defect and study the magnetic field effect on the thermodynamics of the model 
on the basis of the NRG method. Here, as a first step, we consider only the Zeeman field and 
neglect the orbital magnetic effect on the conduction electrons. 



2. Model Hamiltonian with the Zeeman field 

The effective Hamiltonian of the graphene with a point defect has been discussed in ref. [3J, 
and the one-dimensional energy representation for the conduction electrons [S |9] which is 
necessary to apply the NRG method has been derived. In this representation, the effective 
Hamiltonian is written as 

H = ilgra + H dei + H hyh , (1) 

where iJ gra , Hd e {, and -ffhyb are Hamiltonians for the conduction electrons, the defect sp 2 orbital, 
and the hybridization, respectively, with 

H gra = V / de{e - fi)4 a c £a , (2) 
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Here, c ecr is an annihilation operator of the conduction electrons with energy e and spin a =t, 
J,, which satisfies the anti-commutation relation |c ecr , c^ /(j , | = 5 aa >5(e — e') [U EE [9] . ^( x ) ^ s 
the Dirac delta function. D is a cut-off energy that is the same order as the band width of the 
conduction electrons. d a is an annihilation operator of an electron on the defect sp 2 orbital with 

spin a =t, I, which satisfies jffo-,^, j = 5 aa > . fi is the chemical potential of electrons. e sp 2, U, 

and h are the energy level of the defect sp 2 orbital measured from that of p z orbital on carbon 
atom, Coulomb interaction, and Zeeman field on the defect sp 2 orbital, respectively. V is the 
amplitude of the hybridization between conduction electrons and the sp 2 orbital. 

This model of eqs. JI])-([H) is known as the single-channel pseudogap Anderson model |10j . 
where the density of states (DOS) of the conduction electron is proportional to |e|. The Kondo 
problem in this pseudogap system was studied in detail for the case with n = and h = [10] . 
The main feature in this case is as follows. In the particle-hole symmetric (2e sp 2 + U = 0) case, 
no Kondo screening occurs, i.e. the localized moment is not screened by conduction electrons. 
On the other hand, in the particle-hole asymmetric (2e sp 2 + U ^ 0) case, the localized moment 
is screened when the coupling constant V is larger than a critical value V c . In this paper, we 
assume the particle- hole asymmetry and study the case with h ^ and both /i = and p, ^ 0. 

3. Numerical renormalization group method 

We apply the NRG method for the model dJJ-([2]), following the formalism of ref. |1U] to 
investigate the electronic state at low temperatures. The NRG method can include the effects 
of the pseudogap DOS, a finite chemical potential, and Zeeman field. In the NRG method, the 
effective Hamiltonian is transformed into the one-dimensional form with a chain length N, 

H= lim olA-W-WHn, (5) 
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Figure 1. (a) Temperature dependence of entropy of the electron on the defect sp 2 orbital, Sdef > 
under the Zeeman field h/D = 10~ 4 at chemical potential \i = for the hybridization V/D = 0.4, 
0.6330, and 0.8. Here, Boltzmann constant ks = 1. Inset shows S^ef at V/D = 0.6330 for 
10~ 12 < T/D < 10 _1 . (b) Temperature dependence of expectation value of S 2 , (S 2 )d e {, for the 
same parameters as in Fig. 1(a). 



Here, c na is an annihilation operator of the conduction electrons on the nth site of the chain. 
A(> 1) is a logarithmic discretization parameter and a = (1 + A _1 )/2. In the following, A is 
set to be A = 2. is the "hopping parameter", and en is the "on-site potential", and F is 
a factor. These parameters reflect the DOS of the conduction electrons, and have numerical 
values of order of unity [TU] . 

By diagonalizing this Hamiltonian iteratively N times and by eliminating higher energy states, 
we obtain the energy eigenvalues and the eigenstates. Here 300 states are retained. By using 
these eigenvalues and eigenstates, thermodynamic quantities are calculated [BJ. In this paper, 
the localized moment contributions to the entropy, Sdef> an d the expectation value of square of 
the z component of localized moment, (5j)def> are calculated. 

4. Numerical results 

In order to clarify the low-temperature states of this model, we calculate the temperature 
dependence of Sdef an d (S 2 )def m the NRG method. The parameters D, e sp 2, and U are assumed 



to be D ~ 8 eV, e sp 2 ~ — 1 eV, and U ~ 3 eV. Thus we use the parameters e sp 2/D = —0.125, and 
U /D = 0.375, in the following. With these parameters, in the case of no Zeeman field, h/D = 0, 
it has been shown [3] that the critical hybridization is V c /D = 0.6311. At fx = 0, the entropy 
of the defect is Sdef = In 2, In 3, and 0, for V < V c , V = V C1 and V > V c , respectively at low 
temperatures. These three cases correspond to the local moment state, the valence-fluctuation 
fixed point, and the Kondo-Yosida singlet state, respectively. Here, the Boltzmann constant is 
set to be ks = 1. 

Figure 1(a) shows S^ef as a function of temperature for the hybridization V/D = 0.4, 0.6330, 
and 0.8 under the Zeeman field of h/D = 10~ 4 at li = 0. For the case with V/D = 0.4 and 
0.8, Sdef — at low temperatures. However, their temperature dependences are different: At 
V/D = 0.4, for T ~ 10 -4 Z)(~ h), S^ef changes from In 2 to In 1 = 0. Thus, it is suggested 
that the ground state is a spin polarized state. At V/D = 0.8, on the other hand, Sdef changes 
from In 3 to In 1 for T/D ~ 10 -3 . This temperature is much larger than the Zeeman field 
(h/D = 10~ 4 ). Thus, the ground state is considered to be the Kondo-Yosida singlet state. At 
a critical point V/D = 0.6330, we find that Sd e f approaches In 2 at the lowest temperature as 
shown in the inset of Fig. 1(a). 

In order to obtain more information on the electronic states at low temperatures, we calculate 
(Sz)def- Figure 1(b) shows (S|)d e f as a function of temperature for the same parameters. For 
the case with V/D = 0.4 and 0.8, we find (Sf)def — 1/4 and 0, respectively at low temperature, 
which confirms the ground states suggested above: the spin polarized state for V/D = 0.4, and 
the Kondo-Yosida singlet state for V/D = 0.8. It is also found that (Sf)dcf suddenly changes 
from 1/4 to at V/D = 0.6330 when we change V. Thus, V/D = 0.6330 ~ V c /D is a critical 
value which separates the spin polarized state and the Kondo-Yosida singlet state. At V — V c , 
(S 2 z ) de{ ~ 1/8. 

The critical value in the Zeeman field, V c , is larger than the critical value V c at h = 0, 
because the Kondo-Yosida singlet state is suppressed by the Zeeman field. We suppose that 
V c is determined by the relation Tk ~ h, where Tk is a Kondo temperature. The values of 
Sdef = ha 2 and (Sf)def = 1/8 at V c imply that there are two degenerate states at this point. We 
expect that these two states are up-spin state and empty (unoccupied) state, since (S|)d e f = 1/8 
can be interpreted as [(1/2) 2 + 0]/2 [11]. The detail of this unusual fixed point will be discussed 
elsewhere. 

We also calculate the finite chemical potential case, /j / 0, under the Zeeman field. Sd e f 
is always zero at low temperatures. At low temperatures, (S 2 )def; shows the change from 
(Sf)def = 1/4 to 0, as V increases. We do not find any unusual fixed point in this case. 

5. Conclusion 

In summary, the effect of the Zeeman field on the defect-induced Kondo effect in graphene has 
been investigated. The effective model with the Zeeman field has been derived and has been 
analyzed on the basis of numerical renormalization group method. It has been found that, at 
/i = 0, there is an unusual fixed point with the entropy of fee In 2 and the expectation value of 
S 2 of 1/8 at the border of the spin polarized state and the Kondo-Yosida singlet state. 
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